Finding low-energy conformations of lattice protein models by quantum annealing 

Alejandro Perdomo-Orti;,* Neil Dickson,^ Marshall Drew-Brook,^ Geordie Ro5e,^ and 
Alan Aspuru-Guzik^'0 

Dtpaj-tmcnt of Chemistjy and Qif^mical Bioiofjy. Harvard Uitit'trsiiy. 
12 Oxford Street, Cambridge. MA 02138, USA 
'^'^D-Waac Sysitniti. Inc.. 100-44^^ -^^'^^ Crffek Drive. Bumaby, 
British. Columbia V5C 6G9, Canada 

(Dated: 3 Mai'ch 2013) 



'Corrcspoiidinii antlicr'a e-niflil: 



aspuru^chpiiiistry.harvBid , edu 



Lattice proteiti folding inudols are^ a corn^rsCone of computatiunal biuphysicJ^ 
Although these models are a coarse grained representation^ they provide use- 
ful insight into the energy landscape of natural protein^^^. Finding low-energy 
three-dimensional structures is an intractable problenP^ even in the siinplost 
model, the Hydrophobic-Polar (HP) model, Exlianstive search of all possible 
global miniina is limited to sequences in the tens of amino acid^. Description 
of protein-like properties are more accurately described by gcneridized models, 
such as the one proposed by Miya^awa and J ernigaiP (M J ) , whith explicitly take 
into account the unique interactions among all 20 amino acids. There is theo- 
reticaP^^ and expcriiriontaP^ evidence of the advantage of solving- classical op- 
tiiJiization probleiJia uaing quantuiti annealing ^ ^ ^ ' ^ ^ ' ^ ^ over its classical anaioguc 
(simnlated aiinealin^^. In tins report, we present a benchmark implement at ion 
of quantmn annealing for a biopliysical problem (six different experiments up 
to 81 superconducting quantum bits). Although the cases presented hero can 
be solved in a classical computer, we present the first implementation of lat- 
tice protein folding on a quantum device under the Miyazawa-Jernigan model. 
Tills paves the way towards studying optLmization problems in biophysics and 
statistical mechanics using quantmn devices. 

The nearch fut mote ofHciciit eptimizatieii algeritliiun in aii iiiipertiiiit t^iidcfL^Tjr with 
pr(?valent(! on mauy diaciphucs ranging £reiu the social Hciciices to th(? physiciil imd Jiatiiral 
sciences. Belonging to the latter, the protein folding pruljIcnP"^ consists of Rudiiig the 
lowest &'ee-[?ncrgy configuration or, equivalontly, the native structure of a protein given its 
amino acid sequunee. Knowing how proteins fold elucidate thfir threc-tlinLejisioual atructrue- 
function relationship wliitli iti cmcifJ to the understanding ol cuxyiueH mid for the treatment 
of miafolded-protein diseases srjdi Jis Alzlitsiiiici's, Huntiugtoii'&i mid Parkiaaon'a disease. 
Due to the high computational cost of modchng proteins in atouiiatie detaiP^^, coartie- 
gtained descriptions of the protein folding problem, such as those formd Iti lattice niotlela, 
provide valuable hisight about the folduig meclianismiP^ l^^. 

HarueHsing quant uui-uiechajiieal efiecta to speed up the solvuig of tlassical optimisation 
problcTiis is at the heart of quantum aunealing algorithms (QAJiiii^ii^il^. lu QA, quantmn 
methanieal tunneling allows for more cfiicieut exploration of difficult potential energj^ land- 
scapes such as that of claiHie-d spin-gl-iss problems. In our implementation of lattice folding, 
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[luiuiluni fluctuatioofi (tmmelingj Gccurn bctwt^uii ittitus rt?ptesc?nliiig diffc^rout mode! protein 
coufoniiiiliLniJi or folds. 

The tlieorcticiil challtaigo is to cificieiitly map the liard computatioEijd problem of uitorc'st 
(e.g., lattice foldiug) to a clasaicaJ apin-glEifia Hamiltoiiian: such mapping requiring a poly- 
nomial number of (juajitiim bits (qubits) xviLli the nixe uf tiio problem (protein length] is 
doHcribed clscwheit^. Hure we prest^ut a new mapping wtiidi, due to ita expuncntial sealing 
vnth problem size, ih not intt^nded for large iiistanees. The proponed mapping empioyw I'ery 
few qnbita fur small problem iuataneea, making it ideal for this first experimental demon- 
stration and implementation on eurr(?nt quantum deviee^P. A combination of the existing 
pulynomial mapping and more advanced quantum devices wuidd allow fur the simiilatiun 
of mneh larger iiistanees of lattice folding and other related optimisation probleiUK. 

SulviiJR -u-bitrary problem iiistajieea reqiiircs a programmable qiiiuitiuu device to im- 
plement the correapondijig elasnical Hamiltuiiijtn. We emplu\' quantum amieahng on the 
prograjmnable device to obtain low-energy cojiformations of the protein model. We CTnpha- 
size that iiuthiiij?^ quantum mechanical is imphed about the proteiji or its folding procesa; 
rather quantum fluctuations are a tool we use to solve the optimization problem. 

The QA protocol pcTformed here in bIbo known as adiabatie quantum eomputation 
{AQCj^^^. Of all the qiiautmn-compiitational models, AQC ia perhaps the most naturally 
suited for studyuig and solving optijnizatiou problem?P*^. For the experimenta presented 
here, the small finite tempcratiut^ of the aupereonductinp; device is enough to make the 
process less coherent tlian the original formidation of AQC, where the theoretical limit of 
zero tcmperatme and quaai-adiabaticity are usually assnniedi^i^. As we show irj the dii^ 
euaaiou, numerical simulatiojis uicluding these unavoidable environmental effects accurately 
reproduce our ejtpcrimcntal results. 

Experimental implementations of QA or AQC are limited either by the numbei of qnbits 
avaflable in state-of-the-art quantiun devieea or by the programrnabihty reqiiircd to fulfill 
the problem specification. For example, the first reahzation of AQC was performed on a 
three-qubit NMR quantmu devicJ^ and newer NMR implemeiitationa involve four qubit 
experiment Ji^. Other experimental realizations of spin systems have been based on mea- 
suring bulk magnetization properties of the aystems in which there is no control over the 
individual spina and the coupfings among then r^ l '^^ l Quantmu architecturea uaing superccn- 
dueting qubits offer promising device scalabifity while maintaining the abihty to control 
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individual qubits and the atrt^uglli of their interaction tciuplingti. During the prt^p-irEitioii of 
this manuscript, an B4-qubil experiiueiilal Jetermiaation of Ramaey aumbers with, qiuuitiuu 
ajmealing was porfonuetP^j uuderat^oring Lln' progriuinnable capabilities uf the dc?vice for 
problems with over 80 qubits. In Lliin letter, we present a quantiuii aunealijig experiineiital 
implementation of lattice protein models with general (Miyazawa-JemigatP) intcractioiis 
ajnong the amino acids. Eveu though the eases presented here still eau be solved on a 
ckaaical eompiiLc'r by exact enumeration (the aix-amino acid problem has only 40 possible 
configurations), it in remarkablt? that the device anjieals to the gromid state of a wearch 
spaee of 2^* possible eomputational outcomes. This study provides a proof-of-priiidple that 
optiiuiaation of biophysical problems sucJi an protein folding can be studied using quautmn 
mechiujical devices. 

The quautimi hardware emplo^'cd eonsinta of 16 imits of a recently characterized cight- 
qubit unit celP^H Poat-fabrieatioii eharacteriaation determined that ouly 115 qubits oiit of 
the 12R qubit array can be rchably used for computation {see Fig.[l]l. The array of coupled 
superconducting flux qubits is. effectively, an artificial Ising apin system with prograimiiable 
spin-spin coupling and transverse magnetic fields. It is denigiuxl to solve instances of the 
following (NP-liarJ^ classical optimization problem: Given a set of local longitudinal fields 
{hi} and an interaction matrix find the asaignmcut s' = tils', - ■ ■ .y'^. that niijiimize^ 

the objective function ^^(s), where^ 

E{s)= ^ J„s,^,, (1) 

l<j<.V ■I^J<J^-iV 

< li 1^1 < 1, aJid>!, e {-Hi.-i}^ 
Fuiding the optimal s' ia equivalent to finding the ground state of the eorreapondiiig Ising 
classical H^uniltonian, 

jV N 

where trf are Paidi matrices acting oji the ith spin. 

Experimentally, the time-dependent quantum Hamihouian implemented In tlie supercondueting- 
qnbit array is given by, 

H{t) = A{r)Ht + B{T)Hp, r = t/U,,, (3) 

'A'ith f^b = responsible for (juantinu tuunehng among the locah^ed classical states, 

which coricspond to the eigcjistates of Hp (the eomputational b-iais). The time-dependent 
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functions A{r) mid B{t) aro sudi that A{Q} :^ B{Q) mid A{1) <^ B{1)- in Fig.[2);b), wti plot 
theae fuuctiona ii» iiiiplcmcrLtod iii Oil' experiment, i^^^ deootea the time clsipjicd bctwcea 
tin? prt^pairiUoii el Uiu iuiliEi] Htfitc! mid tlie meaauremejit. 

QA espiuit.-i the adiabatic th(?orem of quantum medianits, which stfLtc?Pi that a. quantum 
ayatcm iuitialiaed in the Rtouud state of a time-dependent Hajniltuiiiim remains in tlie in- 
stajLtmieous groimd atate, as long as it ia driven sufficiently slowly. Since the ground state of 
Hp encodes the solution to the optimization problem, the idea behijid QA is to jidiabatieally 
prepiue tins ground state by initialising the quautuni system in the easy-te-prepare ground 
state o£ Hf,, which corresponds to a superposition of all 2^ states of the compiitational ba^ia. 
The ayatcm is driven slowly to the problem Haniiltoniau. H{r = 1) ^ Hp. Dcviationa frejn 
the gruimd-state arc expected due to deviations from adiabaticity, as well as thermal noise 
ajid imperfect ions in the implementation of the Hiuniltenian. 

The first challenge of the experimental iinplenientatiu]i is to map the computational 
problem of interent ijito the binary quadratic expression (Eq, [2]), which we outhne next. 
In lattice folding, the scqiienec of amino acids defining the protein is viewed as a scqiiencc 
of beads (amino acids) connected by strings (peptide bunds). Tliia bead chain occupied 
points on a two- or tbrce-dimensional lattice. A vsilid cunfiRuration is a self-avoiding walk 
on the lattice and its energy is calculated from the sum of interaction energies between 
nearest non-bonded neighbors on the lattice. By the thennodynainic hypothesis of protein 
foldhk^, the global mhiimnin of the free-energy function is conjectiu-ed to be the native 
functional cunfurmatioii of the protein. 

The hydropliobic-polar (HP) model ia one of the simplest possible models for lattice 
foldui^. In tliia model, the amino acids are claasified into two groups, hydrophobic (H) 
and polar (P). To deseiibe leal protein energy landscapes a more elaborate description needs 
to be considered, such as the Mijaaawa-Jernigau (MJ) modeP whicli assigns the interaction 
energies for pairwise intci-ictions among all twenty amino acida. The formulatioji we iised 
ia general enough to take into account arbitrarv interaction matrices for lattice models in 
two aiKl tliree dimensio]i^. In p:uticular, we solved a MJ model hi 2D. the six amino-iLcid 
sequence of Proline- Seruie-Valine-Lysuie-Methionine- Alanine (PSVKMA in the orie-letter 
anilno-acid sequence notation). We solved the pioblem under two different experimental 
sdicmes (see Scheinea 2 and 3 in Fig. each roquirhig a different number of resourcea. 
Solving the pioblem in one proposed experiment-d lealization (Scheme 1) requires more 
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rL^hLiiircL's than the uinubur pf qubits iWidlablc (115 qubitti) in tiiu dcvitc. Sclioiiit' 2 and 
3 axe cxcuiipl(?& of Llif divide- and-coiiqucr atmti^gy, m wliicli aiit pcirtitious the problem 
in smaller iiiHtajicoji and tonibia[?H the irid[?poiid[?nt set of results, thereby obtfiiniiig the 
same yohition for the intractable problem, in the SI aettioji, w<? complement thcae four MJ 
related expeiimeuts with two small tctrapcptide instances {effectively HP model inatancea) 
for a total of nix difForent problem HamiltoniaiiH. We us(?d the largest of theae two instances 
{aji a qubit expcrimeul) for direct theoretical aimiilatiou of the aunetLhiig d\Tiajnics of the 
device. The results from our experiinejit ajid tlie theoretical model, whidi does not uae auy 
adjustable parameters (ail are extracted experimcjitally from the device), are in excellent 
agreement (see panel (b), Fig. 2 of the SI material). 

To represent each of the possible A^-amijio-acid cojifigiirations (folds) hi the lattice, we 
encode the direction of each snceeHsive bond bet^^eii amino acids; thus, for e\'ery A'-bead 
sequeucc we need to apecifv^ - 1 turns corrcspoudinj^ Lo t he iimiiber of bouds. Fbr the case 
of a two diineJisioual lattice, a bond can take auy of four possible directions; therefore, two 
bits per bond arc required l,o uniquely determine a direction. More specifically, if a bond 
points upwards, we write "11"^. If it points downwards, leftwards or lightwards, \vu write 
t)0" , '10'' , or "01" respectively. Fixing the direction of the fiist boiid reduces the (leseription 
of any TV-bead fold to = 2{N — 2) binary variables, without loss of generality. Ah nhown 
in Fig.[2|;a), in the absence of external cojistraints other than those imposed by the primary 
anmio acid sequence [see SI for an example with external eonstrauits) , we can tut the third 
Mnaty variable to tl", forcing the tliird amino add to go either straight or downward and 
reducing the number of neetled variables to £ = 2A' — 5. This constraint reduces the solution 
space by removing eoidbnnations wliieh are degenerate due to rotational symmetry. Thus, 
a particulaj fold is imiquely defined ]yv, 

q = -^01^ t^^l 1-2fli ■ ■ ■ fhN-^^2N-^ (^) 
iurrtl r„,.||j JnrFjl turTilA^-l) 

An example of tliis encoding for a six-amino-acid sequence is represented in Fig, |2|[a). 

Using this mapping to translate between the amino ticid chain in the lattice jmd the 
2(A' — 1) string of bits, wii constructed the energy fimctio:i E{q) in which q denotes the 
remaining 2A^ — 5 binaiy ^miables. Additionally, we penalized folds whitli. exhibit two amino 
acids on top of each other, to fa\'or self-avoidiiig walk couBguratiojis. The energy penalty 
chosen for each problem was sufficient to push the energy of invahd folds outside of the 
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ojicr^ raiigL' of valid tonfiguratiooH (thuao with E < 0). Finally, wo took into aot'oiiiit the 
interaction energy among tlic diffcrcjit amino acids. A dctailotl construction of onr energy 
Rinction for the general CiiHC of N auiino acidH with arbitrary' iiitoraetions is given elsewherd^. 

The experiment tontiitita of the following stepH: a) eonstructioji of the energy function 
to be minimized in terms of the tnm encoding; b) reduction of the energy expreaaion to a 
two-body Hfljoiltonian; and finally, c) embedding in the device. Thewc last two atepti need 
additional rcsoiurces aa expl^unetl beluw. We will foeiia on the siiupleat cxjaiiple [ExperiiuenL 
3; Fig. |3[| to jiliow the procedure in detail The embeddinga for the other five cxperiuieuts 
are provided in the SI material. The energy function for ExperiniCTit 3, cautainiiig the 
contributiojis due to on-j^ite penalties for o\'erlappuig amino acids, and jjair^'iae interaetiojis 
between amino acids is, 

E{q} = E^,'^^ = -1- 4^, + + f)q,^;, - l^jq^q^ (5) 

where 91O (^^^a) encodes the orientation of the fttuith (fifth) bond {sec Fig. |3|- From Eq. [3] 
one vcnfy by substitution that the eight possible three-bit^variable aHsigmnents provide 
the desired energy landscape: the six conformationH with £ < shown in blue hi Fig. |3) 

Eq. |5| describes the energ;v^ landscape ttf eoufigurations but it is not qmte ready for the 
device. Expcruneiitallyj we can specify up to two- body spin interactions, a^a^. and therefore, 
we need to convert tliis cubic energy function (Eq. |5]) into a quadratic form resembling Eq. |l| 
(see SI for details). The resulting expression is 

where the original binary variables and spin opeiatois are related by — ^ (1 — af)/2. 
Experimental measurements of of yield a, = +1 (&j = —1) eorresponding to = (9, = 1). 
Since q, = {i — ^rj)/2, measurement of a[, ^ry, find allo^\^ ils to leconstruet the bit stiing 
qiOqili which encodes the desired fold. 

One ancilla variable wa^i added during the transfbrmatioji of the tlu^ee- variable cubic 
Hamiltouian uito this quadiatii four-variable expression. The meaning of the original vari- 
ables fill f^i, and remains the same, allowing for the reco:istruetion of the folJs. The 
energy of this four-variable espressioji will not change as lojig as the measurements of 
througli result iu values for q]<fi<i:iq^ satisfying = q^q^,. Tliis transformation ensures an 
energy penalty whenever this conditioji is violated. 
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Thu iirdiitcotiirc of the chip liitka aufficioiit ct)Ui:uLlivit\ bclwccii Llic jiiipLrtxinductiug 
rings far ^ aiic-ty-uuc aaaigjiuient of variablea ta qubilH {aoo Fig. Tu satiafv the couudt- 
tivity r(?qiijromt?JitH of tlio four-vari-iblG energy functiou. the cuuplings of ojie of Lho Jiiost 
(.■ouiit^ctod variii.bl(?K| ^^j, wore fulfillod by dupUcating tluH variable iofiide the devici? sudi that 
q.i qi and qi'. In the funn of Et|. [2| the final expresaioii rcprcacTiting the energy functiou 
of Experiment 3 is j-ivt^n by, 

Hp = (Ttjlof -f- + Sff^ + Qcrfofi + 9ff^cr.^, - 2cr| - IGaftrl 

This expression aalisfiea all reciuirenients for tlic problem Haniiltoiiiiui {Ef\. the complo 
tion of wliieh allows for the meaaiirement of the eneigetie minimum eouformatiou of this 
snuJI peptide histanee, The embedding of Eq. |7| into the liiudwaj-e ia shown hi Fig. |4| where 
we label the five qubits used, qi, q-j, q-i, q^, and^j'. Since we want the two qubita reprericnting 

Ui end up with the same value, we apply the maximmn ferromagnetic ctjuplijig (,/ = — 1) 
between them, wliieh adds a penalty whenever this equality is violated (last term in Etj. |7|l. 
These maximmu couplings are indicated in Fig. [T| by heavy lines. The tliinner hnes tihow 
the remaining eoupliiigs used to reaHae the quadratic terms in Eq. [7) color coded according 
to the sign of the interaction ajid ita thickness representing their strength. Note tliat e^-ery 
quadratic terjn in Eq.[7|has a corresponding coiLpler. Hereafter, wc will denote the outcome 
of the flve-qubit me^isurementa as qcspo = OlOOiO^iOg^'/jsl^.ii^.iJ; with = Q {q, = I] whenever 
H, = 1 (Hj = —1). Notice that only the bits preccdijig the divider diaracter | contain physical 
information. These are the ones shown undei each of the protein fold drawing associated 
'Aith Experuuent 3 (see Fig. 

Similar embeddiiig procedurea to the one previously described were used for tlie larger ex- 
periments. Fbr example, in Experiment 1, only 5 qubits define solutions of the computational 
problem. We needed 5 auxifiary qubits to transfonn the c^presaion with 5-body interac- 
tions into an expreji&ion with only 2-body interactions. Embedding of tliia final expression 
required an additional of IS tjubita to satiafy the liardwiu^e connectivity requirements, for a 
total of 2S qubits. Table I in the SI mateiial ammnaiiiea Llie number of qubits required in 
each Btep through to the final ex[>erimental realiaationa. 

Even though the quantum device follows a quantum armeaUng protocol, the odda of 
measuring the ground atate arc not necessarily high. For example, in the SI qiibit experi- 
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juuiit. oiilv 13 uiit of 10,000 iJiCiiBureuiefita yielded t]io dt^ired solution. Wc attribute these 
luw-pereeiLtage-^ lo thu ^alog nature of the device and Lo prctisiou limitations in the rcid 
values of the local fields and couplinp among the qubils in the experimental setup. When 
compared to other problem implementationti, physical problems sneh aa lattice folding lack 
the structure of the Ramacv number ptoblcrrPI In the lattice folding problem implemented 
here, the parameters defining the problem instances are arbitrary and Jo not fall into certain 
integral distiTiet valuer us in the c-ihc of the RaiiiHCy number expeiiment. making precision 
iBSues more prononnced in our implementation. 

To gain insights into the dynamics and e\T)lution of the quantum sj-stem, we nmnerically 
simulated the superconducting array with a Bloeli-Rcdhcld model of the 8-qubit experiment 
(see SI material) which takes into account thermal fluctuations in the states due to the finite 
temperatiiro (20mK] of the quantmn device. For this S-qubit experiment, the simulation 
predicted a ground state probability of 80.7 %, in excellent agreement with the experimen- 
tally observed value (80.3%). It ia important to note that no adjustable parameters were 
used in our simulatiojis to fit the data and all the parametera correspond to values measmed 
directly from the quantum deviea. More details about the iiimierical simulations can be 
fbmid ui tile SI, 



As seen in Fig. |2(e], the temperature of the device ia comparable with the inuiimmn gap 
o£ the eight-qubit Hamiltonian. Therefore, we expect stronger excitation /relaxation near 
the gap cloaiiigi t ^ 0.6, due to exchange of energy with the environment, when compared 
to the other regimen of the annealing sdiedule where the gap is nmch larger than kuT. In 
the abaeiice of environment [a fidly coherent process), our simulations indicate that that 
the Huceeas probability woidd be 1009^1 witliiu numerical error. Fig. |2[d) shows that for the 
simulations at 20mK, the probablflty in the ground state goes down to 55^., but the same 
fluctuations make the system relax back to the ground state, yielding tau 80.27^ suceeas 
probabihty. This is due to the advantageous natural tendency of the system to approach a 
thermal equflibriiun which favors the ground state after crossing the mininmm energy gap. 
As previously (lisciisaed in similar numeiical simulations of quantum annealing algorithm^, 
strojig eoupliiig to the bath ajid non-Markovianity viDuld require going beyond the Bloeii- 
Rcdfield modcP^, but the agreement between experimental and simulated results support 
the validity of the quantmn medianical model uaed to describe the device. Previously 
repotted tcmperatinre dependence predictions for the tuimcling rate on the same qubitsl^ |3] 
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EUld tjxcollt^ut ugrcL^iiiL'iiI, with the sauio luvol uf thecity iis(?d Iiuil' rt^iiilbrto the: validity of oiii 
simulations for this S-qubit iiistaiift-H. 

Wo prL'S(?nt the first qiiajitiuu-much-miL'id iniiiloiuuiitii.tioii of latlico proLoiii inuik'ls iisiii^ 
El prugrEUumablc quantuiu duvicc. Wt? wert? able to oticudt; aud to solve thi3 global uiiuiiuji 
solution for a suiciU tt^trapoptidc Eind hcxapcptidc chain uiidci HCveial c■J^pc^i^lelltal aehGmes 
involving 5 and S qiibits for the four- amino- acid ncqiicnct; fHydiophobic-PoliU^ jnodt^l) find 
5, 27, 2S, and SI qiibiLt cxpciiniciit!^ for six aniino-acid sequence nndcr Llie MiyazawEi- 
Jemigan model for general pairwist: interactioiia. For the experiment wiUi S qubits. 
simulated the dynamics of tht? quiintiim device with a Rcdficld equation with no adjustable 
parfimeterH, obtaining excellent agreement with experiment. Sinte the quantum annealing 
alj-oritlun not oiJy finds the ground state but also the low-lying excited states, it provides 
infonuation about the relevant minimum energy compact structures of protehi tiequeneefP^ 
and it is uscfid to evaluate designability and atability sueb -is that fomid in natiiral protein 
sequences, where tbe global minimum of free energy is well separated in energy from other 
jnisfolded titatctP^. The approach employed here can be extended to treat othei pioblems 
in Ijiophysics and statisticid mechanics such as molecular recognition, protein design, and 
sequence aligmnenP^. 
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FIG. 1. Tiie array of sup £?rcoiid\if ting qiiantimi bits Ie arranged in 4 X 4 unit cells that conEiat of 8 
quantum hits each. Within a unit cell, each of the 4 qubits in the left-hand partitiau (LHP) counectE 
to all 4 qubitE m the right-hand partition (fiHP), and vice versa. A qubit in the LHP (flHP) also 
connectE to the correEponding qubit in the LHP [RHP} of the units cells above and below {to the left 
and right of) it. (a) Qubits ai'e labeled fcoui to 127 and edges between qubits represent couplets 
with, prograrrnnable couphng strengtlis. Grey qnbits indicate the 115 usable qubits. while \'acancies 
indicate qnbits under caUbratioii which were not used. The larger experiments [Experiments 1,2, 
and 4) were performed on this chip, while the three remaining smaller experiments were rim on 
other chips with the same ardiitecture. (b) Embedding and qubit coimectivity for Experiment 
4, coloring the SI qubits used in the experiment. Nodes with the .same eolor represent the same 
logical qubit from the original 19-qubit Ising-like Hamiltonian resulting from the energj' function 
associated with Experiment 4 (see SI material for details). This embedding alms to fulfill the 
arbitrary connectivity of the Ising expression and allows for the eonpHug of qubits that are not 
directly coupled in hardware. 
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FIG. 2. (a) Step-by-step construction of the binary representation of lattice protein. Two qubits 
per bond are needed and the bond directions are denoted as "00" (downwards)^ "01" (rightwards)^ 
"10" (leftwards), and "11" (upwards). The example shows one of the possible folds of an arbitrary 
six-amino-acid sequence. Any possible JV-ainino-acid fold can bo represented by a string of vari- 
ables OlOgiq-ig-.i ■ ' ' with f = 2iV — 5. [bjTimc-dependence of the A{t) and B(r) fimctions. 
where r = ^/'ruFj with ^rurj = 148pi. (t) time-dependent spectrum obtained through numerical 
diagonalization, and (d) Bloch-Redfield simulations showhig the time-dependent population in the 
first eight instantaneous eigenstates of the experimentally implemented S-qubil Hamiltonian (Eq.[3j 
with Hp ftttm Eq. 18 in the SI material. In panel (c). for each instantaneous eigenenergy curve vje 
have subtracted the energy of the ground state, eflectively plotting the gap of the seven-lowest- 
excited, states with respect to the ground state (represented by the baseline at zero-energj'). As a 
reference, we show the energy with the device teniperatm'e, which ia compaiable to the minimum 
gap bi?tween the ground and first excited state. In panel (d). populations are ordered in energy 
from top (ground state) to bottom. AltbougL t = t/t 

fTiTi runs from Q to 1. we show the region 
where most of the population changes occur. As expected, this is in the proximity of the ininiTmini 
gap between the ground and first excited state aioimd r 0-4 [see panel(c)|. 
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FIG. 3. (a) Representation of the six-amino acid sequence^ Proliae-SGrine-Valine'Ly^ne' 
Methiomno-Alamne with its respective one-letter sequence notation^ PSVKMA. We use the pair- 
wise neareEt-neiglil)or ^!iyazawa-,Terlligan interaction euergies reported iu Table 3 of Ref. [5| (b) 
Divide and conquer approach sliowing three different .schemes which independently solve the six- 
amino acid sequence PSVKMA on a two-dimensional lattice. We solved the problem under Scheme2 
and 3 (Experiments 1 through 4). ft) Energ:^' landscape for the valid conformations of the PSVKMA 
sequence. Results of the experimentaily-measmed probability outcomes are given as color-coded 
percentages according to each of the experimental reahzations described in panel (b). Percentages 
for states with energy gi^eater than zero are 32.70^0, 59.88%, S.OO^d, and 95.97% for Experiments 
1 thiough 4, respectively. 
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FIG. 4. Graph representatLons of (a) the foiu'-qubit uiienibcidded cnevgj' function (Eq. [gJ and (b) 
the five-qubit expression [Eq. ^ bs was embedded into the quantum hardware. Li graphs (a) 
and fb), each node denotes a qubit and the color and extent of its glow denotes the sign and 
strength of itJi corresponding longitndinal held, !s,. The t^dges represent the iuterflction coupHngs, 
r/jj- where color uidicatcs sign and thickness indicates magnitude. Smce we want the two qubits 
representing qi {q^ and q^') to end up with the same value, we apply the maximum ferromagnetic 
coupling {-J = —1) between them, which adds a penalty whenever this equality is violated. These 
maximum couplings are indicated in the figure by heav>' lines. For the case of Experiment 3, 
the reconstruction of the binary bit stirigs representing the folds in Fig. |3| hom the five-quibt 
experimental measmemeuts can be recovered by qexp:,i — ^^^^^li^'i^l'^l^lAli' ^ with q, = = 1) 
whenever = 1 (n, = —I). 
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Supplement aiy hilbrmatioii 



SUMMARY 

The paper Finding law-enejyy conformations of lattice proteut inodds by quanium- 
anneaiing prcsc?utH the first CKpcruocntaJ aiid Ifirgo&t qiiiiiitiiiii iuuioaliiig cxpcriiiicjit rolatt^d 
to ail aptimi^atiuii problem in the phyaieal seicneca. In Set. |T| we aiiiuiiiiiriae the eonstnic- 
tion of a more succinct version of tlio energy function describing the [?uergj' landscape of 
the six expc>r:iiiental realiaatiouH of the gc^ueraliKed lattices-folding model using Miyiizawa- 
Jcmigfln pairwisG iaiteraetiojis. In Sec. |Tr[ we present the neccaaary steps to trauafonn the 
energy fmietion into an expression which can be readily implemented in the quantum device, 
hi Sec. |TII[ we describe the quajitmu device iised for om experiments Jiud iii Sec. [TV 



we 



give details about the quantum simulations and results used to support the experimental 
outcomes. 



I. TRANSFORMATION OF THE ENERGY FUNCTION OF THE 
LATTICE-FOLDING MODEL INTO THE EXPERIMENTALLY 
REALIZABLE SPIN-GLASS HAMlLTONlAN 

The energy function for the lattice model cab. be obtained as a sum of different contribu- 
tions. 

Ep(q) = E^^,t.{q) + E^^q) + £^,(5) (Si) 

where E^naiSE{<i) penalities configuratlona with overlaps amou^ any two amino acids, Ejn,[q) 
accounts for nearest-neighbor pairwise-interaction energies among non-bo Jided amino acids, 
and Eg^{q) refers to any external potentials other than the ones coming from iuteraetioiis 
among the amino acids defining the protein. For amino acid sequences in vacuo, only Ea,unte 
and Ep^ are needed. The construction of these tlirec-types of energy functions, in 2D and 
in 3D. fi>i ,in iirhitrary nmnber of amino acids ajid uiteractions among them is explained in 
detail in Ref. |I| Hereforth, we will only focus on the case of eneigy fmictions in 2D. 

1 



A. Cusff of thu Bix-amino acid sequonce PSVKMA (Expotiioetits 1-4) 



Fbr convenience, wc reproduce Fig. 3 of the maiu text as Fig. |ST) wLidi illiiatralcs iiud 
defines the six amino- acid aLtiuonto PSVKMA. 

As explained iji the main test, the description of all possible 2D A^-amino-acid fold in 
vaxiuo taji be described by bit string of leiigUi 2(A' - 1), with the first throe bits held 
constant leaving I = 2N — 5 binary variables as the computational variables of the problem, 

turnl ,^„,'2 ^„^n;[ inrn^N-L) 

For the case of = 6 {seqiiejico PSVKMA). the problem Is completely spocifit^d by the bit 
string 

Qi^AA = ,^01^ ^r/i q-ii!:^ iiiijr. ifi.qj . (S3) 
lurnl tjirnS J^irna lTim4 tiirT[5 

Ey using the construction in Rcf. (Tj tlie 7-bit energy function describing the sequence 
PSVKt/lA [Scheme 1 ui Fig. |SB is given by, 



^PSVKMA(5i'A.4) = -ti-i + S^i^i -f- i^qzffi - IR'/i^a^a - Sgiq., + 12^19-^9,1 + 4939^ -f dtjiq^qi 

- ^'hfi-iff^ - 12^i9:^</.|</r. - ir)9;i</.|y=i -f 28yjy;iy,9^^ - 2^,9,9,- - 19;h-/,- -h 2</jqj^Q 
+ IS^iq-jf/jii/r, - 2(/i9i9B -f- 4gi9oy.i(7rj + 2q;[^iyt + 13^i9;i(/i(/^ -f 1(/L'y;igi96 

- 335i^i^;iy5VQ - ^Oqiys^f, H- 13y|y,^|ig^| + Ig-j^j^a^b - S^gi^^^jq^Vc + J^Vi^J^a^O 

+ 2^^jq.,</7 -h 13</iyjy.,q7 + 4q;Hy.i^7 + i^f/j^i^i^ - 33^iq.yj5yiy7 -f- 

- lS?,j^r>qr + ^q-i^hfif-qi - 23q\iijqiq^il7 — 33q±q:iil.iq^q7 -f- 62^i^jg;iy.iy-i[/7 + Tqii^r 
-f- 2q3qo^ + Idqiq^qaqr + -^qaqofl; + 9q:!yi<;{></7 - 33^i<;'j^;!qf,qr - 20^.iQrf,q7 

+ I3^iq.i</fiy7 + ^^Liy i^fl^ - ^Tqiq-^q^iiiaq'i + '^q\iq\qv,q- - ^^qiqi^q iqi^qr - '^'^q-q.^q-v'k^- 
-f- 90^1 1^1,9- - iSq=iqh?r + ^q-iq^q^q-i - '^'^qiqiq^.q^^i - 3^'lq.\q^q'^q^ 

-f- 62yi<;'jy;!qr,qitg7 + fiZq^q^^q^qi - Z'^qiq iq^qsqy - ^Iq-iq-xq^q^q- + ^^qi,qi.q-iqLqiiq- 

- ^^q:iqiq5qi,ii7 + G2qiqiqiq^q,,ii7 + ^^'l-2q:iqiq5qi,q7 - ^^^^qiq^q^iqiqi^q^qY 

(S4) 
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FIG. Si. fa) Representation of the six-anuno acid fieqiience. Proline-SGrine-Valine'Ly^ne' 
Methionine-Alanine with its respective one-letter sequence notation^ PSVKMA. We use the pair- 
wise nearest-neighhor ^!iyazawa-,Tprlligan interaction energies reported in Table 3 of Ref. [5| (b) 
Divide and conquer approach sliowing three different .schemes which independently solve the six- 
amino acid sequence PSVKMA on a two-dimensional lattice. We solved the problem under Scheme2 
and 3 (Experiments 1 through 4). ft) Energy' landscape for the valid conformations of the PSVKMA 
sequence. Results of the experimentally-measmed probability outcomes are given as color-coded 
percentages according to each of the experimental realizations described in panel (b). Percentages 
for states with energy gi^eater than zero are 32.70^0, 59.88%, S.OO^d, and 95.97% for Experiments 
1 thiough 4, respectively. 

As tihowii in Fi^. |ST| oxproasioiiH for c?aL'h of the diflcrout t^xpcriiiiciili in SckciucH 2 jiud 3 
can be scqucntiaJly obtEiincd by fixing the value of some of the variables in ^'pgyKMAi^'i^i^)' 

The enetgy functioji for Elxpcrinieiit 1 Ls ubtiiiiied by^ evjiluatiuj> £psvKMA(^r>.4.4) 

3 



^1 = 1 (lliird amiao-tH:id iiiovt^ Id Uio riE;lil) mid q-, = (foiirUi iiiiiiiKj-iU-iil iiiovn;?! oi- 

tlu;[ (km'ii in rij^ht. oxpl(ji1iii]> iippi;[/Lnvri liHlf-plium; Hvmmotry). AFtor irhibi'liii^ Ike 

fivn; iciii-iiiiiiig viixiiitjlca that thi;ir Librls Tnnii i-'i ifiisU'-ul <il i.r y^, u = 

OlOf/i i/L.j/'^';['/-,r/i,y7 — - > = OU)l()i/]i/:jj/<'|[f/-,; the rcaultuig cxpicasiuu ilcitriLmg 

rfl.il.i I 

the ciicigy liiudscapc fui Expi;i iiiK;iil. 1 ih givi;ii by 

+ 7^4^^ + 4515455 - 7595496 - 2451^2^^55 - l^^q^q^ + 2053535455 

+ 29^1 f^^yj! ^456 

The energy function for Ebqwriment 4 is obtamed by evaluating ^psvKTvtA (^^aa) ^ith 
^1 = (third amino-acid movo?i down). After ri'iiamiiig the hix rt;riiaiiiiiig viLriiiljlo!^ so 
that lhi;ir labeb^ jipiLii 1-6 iiiHtusad of 2-7. i.v.. <1{.a\ = ^^^iiiii2ff:iff^*!^iii^il7 ■ > Qrzp'i = 
U [[)[)/) iij2t}\iq.i<j^q&, the resulting exprej^sioii deiCriLiug the energy landscape for Experiuient 4 
ih gi\'i;n bv 

^psvkma(^^p*) = "^1 + 1^?!^ + - ^iq^q^ + Miq* - ^^qig^q* + i55a5i - ^15^54 

— 155a5354 + 2S5ig25^5j - 45355 + 25153^^1 + 2525^55 + 451535355 + 7545£ 

+ 75556 + 25i54«fi + 4«9g-ig& + Sgi^^g-j^B - 2053545^ + 4515^5^55 + 959?3g45B 

— ZTqunqiif^q^ - ■i^i^f, -\- Ufii!2ih + 7ffiyi> -h 25i5a5fi +45253511 + ^qi'f^q:^^,. 

+ ^Qiq^qm + ^qiq^q^q^ + 95a9ay&^ - ^'^^mq^iqm - iS545&^e + ^qiq^q^q^ 

— 335i5954g&56 + S3^3gjgB5a — 3^5153545556 — 33525^5^5558 + 9951525^5^5559 

(S6) 

The energy function for Experiment 2 is obtained by evaluating ^p^vkma'^^p*^ "^^^ 
Qi = d [fourth amino-acid moves either down or right). After renaming the five remaining 
variables so that their labels span. 1-5 instead of 2-6^ i.e., ^^1714 = 0100515^53545556 ^'""> 
Qexja = 010005153535455, the resulting expression describing the energy landscape for Exper- 
iment 2 is given by 

- 2[}ii2q:iif^\ -h dqiq^q^qA "H l<^iq^ + ^\qiq^ - iSya^5& + 75*^5 + '^lyi^S (^7) 

- 205254^6 + ^qiqiq^q^ - ISgagj^e + 53525^54?^ - 335i5a5a^4«5- 
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FiiirLlly, tliu iL*iiuig>' luiictiuu lur Exputiuioiit 3 iti Gbtaiin^ci by t^valuatuig ■^j^i^Yi^]j\^(^f'-Tpi) 
■^illi (f[ = I. q2 = 0{Foiiith iiuiinu-aciti jnovcs left) iunl tfi = Q {fifth ajniiio-iLCiil luttvc* citlitr 
down or itih), (?Kp]oiti:ig tin- coiiHtraiiis uopOHcd b>' the tlirt'C fixed ajniiio-iicidti [P,S, iiiid V). 
After rcnajuljig the throe rcrnaiuiiig vHii-iblcs so th;it their labels lire yi, ff2 -md inatead 
of g2,q5. inid ^7, i t-'-, q^^pi = OlOO*jiiiiq:iq.ig^.% ^ - '"'"^'""^^"A q^pa = OinOlOyiO^q:!. the 

tcIeldg] 

rcsultiug cxprts&iun deacrihing the energy laudaeapc for Experiuieut 3 ia given by 



B. Case of the foui'-aiiiiiio acid sequence HPPH {Experiiiioiit 5) 

Besides tlie nix-aiiiiiiu iu.id setjueuee euusidered atjove, \ve alao euiiatnitled [.he eiie[g,y 
foiictioii fur the simplest of all sequences within lattice proteiu modelH, the HPPH lour- 
amlno acitl neqiieiice witluii the HP model. For = 4, we ean specify any of its folds by the 
bit string qe^p^ = OlOqiq^iii. The three-bit eiierg>' function describing the energy landscape 
of Experiment 5 [see Fig. is given by, 



C. Case of the four-amioo acid sequence HPPH under external constraints 
(Experiment 6) 

A mure realistic in iiivo picture involves the presence of chaperone pioteina uasisting the 
folding dynamics towards the global minima. Chaperones, molecular docldng, ajid molecular 
recognition are examples of problems which can be studied by atlding external potentials, 
Eext{fi}, beyond the intrinsic interactions defined by the amino-atid chain, Eoji^,ii[q) and 



Epu.{<l) [see Eq. Si l The first consequence of adding an external potential Errfiq) the 



chaperojie-like environment surrounding the small four-amino-ricid sequenee HPPH, iWuh- 



trated in Fig. 32 by the pink-shaded area near the peptide) is that we can no longer exploit 
the s>inmetry of the solution space for upper and lower half plane confonnations. There- 
fora, we cannot set the first variable of the turn associated with the third amino-aeid to 
aero. Under external potentials, we specify arbitrary folds of the foiu:-amino acid problem 
by qexp^ = Olf/iqii^iijii where giqiiQ^tti) encodes the orientation of the second (third) bond. 
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FIG. S2. Energy landscape lor the four ainino-acid seqiieuce HPPH, (a) in vacuo {Expenrrieat 
5), and (b) under the presence of a chaperoue-like envirounieut [Experiment G) reprisented by 
the I'ed-shaded rcyon. lii panel (a) [panel {b)]. percentages indicate the experimentally measiiri^d 
probabihties of each state collected over 10,000 [28,G72j runs of the quautimi annealing algoritlmi 
described in Sec. |lil C\ In the case of Experiment G. numerical results from the Bloch-Redfield model 
discussed in Sec [TVj are included in parenthesis. Although the vaiiables involved in Experiment 
5 (Experiment 6) axe described by qapo = OlO^iflaffil^Jfl^i' (^expfi = 01fli^Jfl;ifl.j|fliiflri3aJ5d')i under 
each fold wc write only the physically-relevant varia,bles which define the conformation. Since we 
show some experimeutal outcomes also for states with E > Q, then it is natmal to find states whicli 
violate either the and condition or the ferromagnetic condition; for these cases we explicitly write 
the aijxilhary ^■ariablES which went into the quantum hardware. For example^ in Experiment 5 
[panel (a)], the state Qetpq = 01fH01|ll violates the AND condition since q.i ^ In the case 

of the state q^rpii = tll0010|01 the ferromagnetic condition for qi is violated since ?t q^-. Each 
overlap of the amino acids with the cbaperone raises energy by four unitSn whereas overlaps (red 
crossings) among amino acids in the choin raise energy by two units. 

The L'XtLTual potouti-Ll pt?JiEi]izc?H tuiifutiLiatiGiiH ill which cither the tliird or foiirtli ainiiio 
add. go into the chjipcrtmc region is: 

E,,,,{q^^) = Xt,{l - - q.,) ^ XUi - q,)^., + X^^^^^ ^^^^^ 

Tho penalty A^^^ raises energy only when tho third amino-acid moves down (qi = 0,g? = 0], 
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Apjj irLincw L'TitTj^y uuly wiluii Lliu tliird iuuuio-atid iiicjvch right {fjj = O.ij-^ = 1), A^^j raihL'H 
cucrgy uiily when thu Uiird aiiiiiio-acid movea dcmii iuid the fourth- cmiiiio acid moves right 
{^1 =0.(j-2 = 0.^;, = fl, = 1), cUid the lEiat penalty, A^^. raiaos eiiorgj' o:i]y when the thhd 
ajnuio-aeid moves dowji and the fourth-aiuuio acid moves right (91 = = l.(f^i = fl, fj^ = 
0). Each overlap of the amino acJds with the chaperoiie increases enerj^y by four nu'iln. i.e., 

\d _ Vf- _ Jifi- \rjf _ I 

When the third amino acid is altio ahowed to mvtve upwards, the energy hinttion for the 
HPPH ehain in vacuo is given by, 

(Sll) 



After adding Etj. |S1Q| aud Eq. |S11[ the residting energy function for the HPPH peptitle in 



the preaeiice of the "^chaporoiie"' envirojunent ihuatrated in Fig. [S2) is given by, 
EifPPiixhipigeTpb) = 1 - 391 + - 4^1^ - 9;h + qiq^ - 2q^qi + Aqi - 2qiq.i 

{Si 2) 

- S^j^i + ^jqiqiqi - 2q:iq^ -f- ^q-^qi^qi - q\q'iq:^q.\, 

n. EMBBEDDING OF PROBLEM INSTANCES INTO THE QUANTUM 
HARDWARE 

A, Reduction of hi)*h-ordcr terms to a 2-body laing-like Haioiltonian 



As cxplahkcd in Liie injdn text. cJthough the above energy expressions (Eqa. S 1 1. ^5 SG ■i. 



S7 ii. J^R ■j, and |Sl2p ) describe the desired energy landscape, they are not yuitable ibr 
experimental implement atioii. We need to reduce the degree of the high-order terms {cubic, 
cuartic, etc) to a quadratic expresBion (up to 2-body iuteractiona). These iiigh-order terms 
indicate m-iJiy-body inteiactions which are not experimentally fea^ble within the cmrent 
quantinn device. To achieve this 'Aithout alteruig the low-energy spectra [E < 0) where the 
target minima is iriupjjosed to be fonml, we use the teelmique described in Ref. [TJ^ . In the 
maui text, we i>reaented the simplest eatie where oidy one reduction was required (expresiion 
for Experiment 3. Eq. In the following we will focus on the next moat complex case 



(Experiment fi. Eq. ST2 1 wliich can be eaaily generalised to obtain any of the 2-body energy 



expressions for the larger experiments. 



We introduto two iuiciliii biiuirv vaiiabloH, i/^ and ^q, and subntilutu Eq. 1^^ 12 1 with products 



of the form ^ig^ mtu aud q-iq^ into g^. Tliia siibatitutioii traiisfuniiH tlic cnaigy expression 
(Eq. |H12p into a qnadratie oxpreasioji, e.g. the liighL'St-ordcr term which is quadratic, qiii-^fiifi^, 
is replaced by t/oqa whidi beeomes quiulratic, as desired. Under these jiubstitiitiuuH, whenever 
we h-Lve six-vflTiable aasigtiments, q^pa = Olgiga^a^jlga^ii- auch that = qi f\ = quii 
and q,-- = q:i A (ji = q:iqi, we ha\^ the same energy spectrum as the one for the original 
quartie, four-variable cxprcjiKiun. Since tlie^e two -uieilla are new vJiriabk;H \vhosc values 
are independent of the four origuial variables, we need to penahae six-variable as^iignments 
whenever ^ qiqi and =^ q^q^. For every "eollapse'" of the foim q,qj — > ri^, we adtl the 
penalty EA{gi.qj, = ^jj(3rit + giqj — 2q,rk — 2qjS,j). where J.j is a positive number 

representing a penalty ehoacji (for more details see Rcf. m auch that ansi^muenta violating 
this AND condition correspond to free-energies E > outside the relevant search region 
{E < 0). The function E^Xq„qj,r^-S^^) = only if = q,q^ arid E^(q,,q^. r^;^^^) > 
If fh 7^ lilj- The fcix-variable expression reaultiug from the insertion of the new aneilla 
variables plus the penalty fmiction bccomea, 

+ EA{qi,q^^q^;Sii) + ^A(9a, ^J- ga^ 

= 4 - 3gj + ^2 + &qiq2 - (f:i + (liq.i - ^q-jq-.i + -Iqi - 2q\qi - ^q^qi -f- Iq.iqi 
+ 1^5 - l^qj^B - 12q;jf/B -f- ^q^q^ -h IQqo + 5q-2% - ^^m. " " 'i.'.^'i 

(313) 

where acconling to the criteria in Rcf fH we have chosen = 6. and ^;!.| = l. 

To rewrite this quadratic form in terms of the spin variables [si], we apply the transfor- 
mation qi = j[l — Si) to each of the binary variables, 

13 3 7 1 3 11 

- 2-^'-i'^i - ^.'^-iJi.i -f- a-^Si — 3*1*5 - 3*^*5 H- -.y.j^fs H- -.y^^s - 2i;^^|^ (Sl-l) 

1 
4 

After Bubstracting the constant (independent term), we can fulfill the requirement that 
< 1 ajid \J^j\ < 1 by sealing all coefficients of Eq. ISl ll du^Ti by the maximimi absolute 



valut? of aJ] cucfBtiyuth. Tkc rt^uunuiJized quadratic CKprcsiiioii is given by, 



= (13.^1 + H- T^a + *4 - — Sss + SsiSa + aisn - 2,y^«a . 
— 2.S1.S., — Sft^^i.i -f- -l&;[ft.i - 12,S|ftr, - 12^2^^ + S.y.i^Ta H- Bft^^ig — S&:iJ^g 



The fiual laing apiii-ghis^ Hjaiiiltuiiian (before cuibcdding into the quantum dciicc) caji be 
obtahiod by the substitutiuu s, — > a^. 



B. Einbc^ddiiig tntu the? quantum biirdwfu-o 

Eq. IS 161 dues not fulfill Uic ciiiij-couuectivity roquircmentH (see Fig. 1^ for the primal 



graph roptcjieutiiig Eq. pTHl This limitation is trxcd a,t the cost of adding two new qubita 



serving as replicaa of the two qubita wliidi arc linked by more than four conu(?et:oiis. To 
enforce that the repHcflfl of the T-th qubit [of.) produce the same outcome as the original 
t-th qubit, we couple of and trf, with a strong ferromagnetic coiipling, such that whenever 
th(? outcomes of the two variables are different they get penalized by a chosen penalty 
factor 7, > 0. The function whicli performs this penaHzation for each leplica f-t\i qubit 
JH i:/F,v(Kr};7,) = 7,(1 - tf^af,). Notice that HFM{[of = 0, if = = ±L but 
HFM({<ti}\ 7t) = 271* il s, ^ a,'. Fbr this atudy, a vjdue of 7^ = 74 = 1 suffices to leave 
assignments wdiidi violate thia condition outside the region of interest with E < 0. 

The redistribution of the coimections ajnong the original ajid primed qubits is given in the 
right p-mcl of Fig. The modified function taking into account the added feiTomagnetic 
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couplings is, 



H- 73 = 1) -H Mp,M({(7=y, 74 = 1) 

= {IZal 4- Su^ + H- - Sa^ - + 6;rfa^, + trja.^ - 2a|tr.^ {S17) 

— 2iTfcrj, — 8£r|,iTj + 4D■^t^,^ — 12tr|£r^ - l^D-^tr^ + ^o\o^ -\- ^o^^a^j - Btf^O^ 

- b^l^^-g - + (1 - a^fT^O + (1 - O-^IO 

Again, wc? subtract tiic; independent constant tirma from the insertion of the H^m functions. 
Tlic final expression, wliieii is iinplcjucntablc in the qnaiitiim device is, 

H,^^5 = -2 = (13(rJ + ^ix^ + 7(7^ + <x^, - ^<fi - ^<xl -f- i^Wi' + ^Wi - H'^Io"! 











Encoding lor Qnpenrnent 5 Anh Ernodlng eKperiment 6 
(HPPH) « w tt « (HPPH wiTh chjpcrone) 

TBfUlTLm! TLiri eAj J*Cfl» -^^^ -yj^ junj 



FfG. S3. In the graphs preEented in (a) and (b), each node denotes a qubit. and the color and 
esctent of its glov denotes the sign and strength of iti eorrespondmg longitudinal field, Aj. The 
edges represent the interaction couplings, where color indicates sign and thickness indicates 
magnitude. The maxiiuuni couplings are indicated in the figure by heavy hnes. (a) Piimal graph 
(left) and the embedded representation of the expression implemented in the quEUJtum hardwaie for 
Experiment 5 (HPPH m vacuo), [b) Primal graph (left) and the embedded eigLt-qnbit expression 



(Eq. yiS ' for E^mment fi (HPPH in tin? chapemue-lke environment). 
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TABLE T. Niunber af qubitEi needed for each one of the six experinieuts described in Fig. ST arid 
Fig.jS^I Thi? moEit compact veiaion of the euevgy tuuttion correspondE to the secoud colimm. Each 
one of the steps, reductiDn of high-order terms in the energy function towards a 2-body Ismg-hke 
Haniiltouian and eiubi^ddin^ of this Isiiig expression to fulhll the physLcai connectivity af the gubits 
in the device, requires more auxiliary qubits. The huai cahiinu reports the number of qnbits in tiic 
ExperLmentailv implemented exptessLon of the energy fimctiou. 





Number of qubits needed 


Experiment # 


energy timet ion 


Isiing HaniiltoniiJii 


hardware-embedded expression 


1 


'j 


10 


28 


1 


.'] 


10 


27 


3 


■i 


4 


Tj 


4 


tl 


19 


SI 


5 


■i 


4 


Tj 


G 


4 
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Tiic cmbcddiuga fur Exptiriiiieiits 3 iuid 4 arc shown in Fig. 4 and Fig. 1 af tlic mixlu 
textj respectively. The embeddhigs corresponding to Exporimcjit 5 and 6 aro represented in 
Fig. |S3] while tiio embedding for the mediiim size problem uistHJices (Expeiinientti I tmd 2) 
Eire represented hi Fig. |S4| 

III, EXPERIMENTAL DETAILS 
A. The processor chip 

All experlme:itH discussed liereiii vL^re conducted on a sample fabricated in a lour Nb layer 
superconducting integrated enemt process employing a standard Nb/AlQx/Nb trilayer, a 
TiPt reaijitor layer, and plauarhied SiOj dielectric layer?; deposited with a plasnia-eidianced 
chemical vapour depositioii pruccs!-i. Desi^-n rides included 0.25 f-iin lines aiitl spaces for 
'Airuig lai'ers and a jniuimuiu jmietio:! diameter of 0.6 ji-ni. Experiments were conducted in 
ail Oxford Intitrunients Triton 400 Cryofree DR at a temperature uf 20 mK. 

The sample processor chip contains a coupled array of 12S qubits uf a (lesign disciissed 
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Experiment 1 



Experiment 2 




FIG. S4. Embedding of RxpGriincnts 1 and 2 into the quantum hardware. The 28 qubLts (27 
qubits] from Exptrinient 1 (Experiment 2) ha.ve been relabeled to show the qubits which were 
strongly ferromagnctically coupled representing the same variable and biased to have the same 
experimental outcomes. Both problem instances resulted in ten-qubit spin-glass Haniiltonians after 
reducing their energj' expressions to the Ising-hke 2-bodv interaction expression. The additional 
qubits are part of the embedding procedure i;sed to fulfill the arbitrary eonneetivity of the Ising 
expression, allowing for eouphngs of qubits that are not directly coupled in hardware. 

iu Rcf. |3i Each qubit is an if-SQUlD flux qubit with duublc-wcU putcjitial, as depicted 
iu Fig. TliiTv' rUC luagOHticidly coupled with sign and uiiiguitudc tuiiablc ouuplorti iu 
a manner dcaciibctl iu Rc£. |3| The array is built up of 16 cight-qubit unit cells. For 
example, Experiment 6 waa conducted using a auigle uuit cell (higlilightcd iu Fig. The 
connectivity of qubits witliin the unit cell is ahowu schcTnatically in Fig. [SHJ^. 

Three dificrent diips available with this ?iame arciiitecture were used Lo nui the different 
problem instances (ExperimentH 1-6, Fig. [^auclj^. Experiments 1. 2, tuid 4 were run iu 
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chip, while ExpcthuLiit 3 and 5 u>pCtl a diflcruiit chip. Exporiiuuiit 6 uncd i.hv naiiiL' eliip 
and unit used in Rcf. [HI Since llic diips li-wi: the same areliitccturc and design but 
different calibration parajaeters, wo will foeiia un the L'liip ijyod to run Experitnejit 6. and 
report all the parameters us(?d to run the uiiiiierioal simulation reported iu Sec. |TV) 
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FIG, S5. (a) lllustratian □£ a fiiu^^le rf-SQDlD flux qubit. is the flux bias appUai to the 
major (qubit) loop, and ^-jj. is thi? flux biaa applied to the minor (CJ.l) loop, (b) Cross-sectiou 
of the double-w£?U potential of an rl'-SQUlD flux qubit, with 4 localhted eui?rgy levels mai'k«i. 
primaiily iiifectii the qubit bias f. whereas aflects both the barrier height SU^ aud t. 



Magnetic Environment 

The maguetie fit^d in the sample space was toutroUed with three coueentric high per- 
mcabihty shields and an iunennost auperconducliug shield, FVuiiher active compensation 
of residua] fields was adiie\T3d with cojupouKatioji coils oriented along tlaac axc&, ajid us(?d 
in conjunction witb on-chip siipetcojiducting quantum intcrferi^uce dei^ce (SQUID) luagrie- 
tomctera loca-tcd Jiear each of the four tomera of the processor block (Fig. Compen- 
sation coils were adjiistcd to minimise the magnetic field measured, at the magnetometers 
w^hile the chip was at -\.2 K. The eliip was then thennafiy cycled just above ajid then bach 
down tliroii^h its superconducting tr-uitiition temperature at this muiimal field. We estimate 
that the chip wiLs cooled Ihiougli its auperconductiiig Lrimi^itiou with a field normal to the 
chip surface \B±\ < 2.5 uT. and that parallel to its surface < 3.6 nT over the area of 
active circuitry 
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FIG, S6- (a) Optical pLoto^rEiph d1' a portion of a pailially fabricated 128 qiibit chip. The block 
of eight qubits used in this i?xpeiimerit is outlined in red, (h) Artificial spins are connected in 
a complete bipartite graph A"4^i, and interact via couplers which are continuously timable from 
ferromagnetic to anti&rroinagnetic interaction. A line between artificial spins indicates that a 
coupler is present. The coloming indicates one possihle anangyment of coupler settings, (c) An 
example of how a linear^ ferromaguetic laing spin chain could be implemented by selectively tuning 
some couplers to a ferromagnetic setting {J < 0) (green), and turning oS the rest ( J = 0)- 

C, ExperiiuGutal luethod 

The experiment discusaMl m the umuiXHoript ia outlined in Tablt? |TT[ 

TABLE 11. Outline of experiment 

L Initialisation 

1. Calibration: measure intrLusic device parameters such as junclion Ic, qubit inductance^ 
transformer mutual inductances, etc. 

2. Homogenisation: use on-chip programmable fiux biases to ensmre Ij, of the diiTerent 
qubits match during annt^aling. 

Tt. Annealing £^ read-out 

1, Set h. J 

2. Anneal (reduce A[t) apH increase B{t)) 
[J. Read state of spins 
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Till' slL'ps ill piirl J wtTu pcrlunLn^d unto iuid would, in ^oin^raJ, ciily bo purforiiicd oiicc 
far ii new tliip. Tlii.' calibration step I-l is performed by measuring the circulating cinrcnt 
!p in oath qiibiL. and its dependence on the CJ.l loop flux biaa Prom tliifc informatioii, 
one can extract the qubit eritieaJ current and ijiductfiJite L, Details of tluH procedure 
discussed in detail in seclicni IV .A of Ref. |31 Given these qubit parameters, the effective 
inter-qubit coupling strength attiiincd by the tunable couplers tan be detcnniued. Tliia was 
done by measuring tlie difference in luagiietit flux coupled into a qubit B between states | 1") 
ajid I l) uf i\ qubir. A. This coupled flux was measmed as a function of the setting ef the 
tunable coupler between qubits A aud B. in a manner described ui detail in Ref. IS 

Once the device parameters for each qubit have been extracted, the effective jimctJon Ic 
Euid inductajite L of each qubit arc tuned with on-chip tuning ntructures bo as to make them 
as siniihu to each other as puasiblc. The goaJ of this humugenisation procedure is to ensure 
that the circulating currents, Ip. of se\^ral qubits remain elowe to each other in niagniLude 
while the qubits imdergo aune-ding. This procedure is discussed in detail in Refs. H and 
On-chip tunuig structures enabhng this homogenisation are also described in Refa, and 
151 Fij^me shown the superimposed plots of the measured circulating ciUT^ent Ip (left) 
ajid tumiel sphtting A{t) (right) of each of the eight qubits uacd in tliia experiment after 
liomogejiisation. Qubit capacitance is extracted by measuring the spacing of macroscopic 
resonant timnelhng rate pCiikiJ^. At any point in the standard deviation of the measured 
Ip across the S qubitK is less than 25 uA. The uncertauity hi each measmi?ment of Ip is about 
9 nA. The homogenised <:le\'ice paTameters are summarised in Table |TTT1 

The steps in. II are perfi^rini.'d re|jeate(lly. Step II- 1 in where the Hamiltonian parameters 
k, and from Eq. ISlfll are progrannned. For each such problem npeciflcation, steps II-2 



aud II-3 were repeated to allow collection of statisticti about the lelarive probabilities of 
the possible states. For data presented in thin paper related to Experiment 6. Il-l xvas 
repeated B times, after each of which^ II-2 and II-3 were repeated 4096 times, for a tot-d of 
32 JfiS repetitioiLri of II-2 and II-3. However, step II- 1 non-negUgibly heated the diip, so in 
order to allow ample time for the chip to cool b-ick to the base temperatiure. the first 512 
repetitions after each execution of step II-l were removed, leaving S x 3^534=28,672 total 
repctitiojis of II-2 and II-3. In the case of Experiments 1-5 the statistics were collected 
ovLT in, nnn measurements hi each experiment and enough thennalization ti:ne was allowed. 
Therefore, -dl data was included, in the statistics without the need for removijig any of the 
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FIG. S7. (Left) Measured circulatiiig current Ip of each of the eight qnbits used in thi& experiment 
after hoinogenlsation [step 1-2). (Right) Comparison of measm'ed tminel .splitting ^(r) (labelled 
Aq) for the eight qubits used in the experiment, and the ^(r) fit to a physical model of the rf- 
SQUID. Fits of measured tumiel splitting j4[r) aje used in conjunction with fits to Jp and MRT 
rate measurements to estimate parameteis shown in Table |ITI| 

iiiltta! moaaur(miciits. The experimental reaulta of the prubabilitiea m[:ii>iur[:d are reported 
Eia percentages in Figs. |BT| and |S2| 

Aimealing waa perfurm(?d by raiyiug the single qubit tunueliiig barrier. Tliin is aoeum- 
pllahed by ehaugiiig 't'-ir linearly in lime, frttni 0.592 Lo 0.652 aver a period of 14S fis, 
Eia aho\Mi in Figure ISSl Cireulating eiu^reiit Jp tJiown in Figiui? ISTI is plotted over exactly this 



range of ^oj.. This also has l,he effect of changing parameters A{t] and B{r) from Eq. (3) of 
the main paper, as yhown in Fig. lb of the main paper, and as discussed in Hef. HfS Control 



points a and ^ in Figiure S8 eonespond to the beginning and ending times of Fig. 2(b) of 
the mtun text. 

After the qubits have completed annealing, when ^llj, has been set to $ii as shown on 
the right in Figure |5S| states of the spins are read wi1h a hysteretie dc-SQUID readout, as 
described in Ref. HH 



D. Tlicriiioiiietry 

In addition to a Rutheniiim Oxide thermometer mounted on the dilution refrigerator 
mixing chamber, the effective qubit device temperature obtained dining the measurements 



TABLE JH. Tbtal Josephson junction critical cnrreut, qnhit inductance, inductance of loop 2, and 
juuctmn capacitance extracted from circulating current and tunnel splitting measurements, and 
Macroscopic Resonant Tunneling (MRTj peak spacing. 

Qubit Li (pH) (pH) C(fF) 
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di&tiiiaaod in llic luaniiaiTipL was dt^tcnuiiicd iii twu i:idopc:idtut WcLvs. The first is baaod uu 
analysis of tho singlE^-qubit Macrostopic Ris^onant Tunnelliug (MRT) rato, and ita dopt?3i- 
dence on. the qubit loop flux bi;is Oi^. MeaaiurcmGiita and analyaia of MRT rates for the 
devices used iu this ciX[>eriiiioiit are discussed in Ref The aeouud is based on lueasureiut^ut 
of the equiHbriuin P\ vs. ^ij. attained at fixed barrier height (fixed value of ^>ir)- Both of 
tli(?ye l,eohiiique?i are dlseussed in some detail in Ref ^ 

At a lixed barrier height achieved with a fixed \'ahic of 'l^i^. the equifibrium probability 

apptoatbea the thennal distribution: 



2 \k 

where Ij, ia the vjilne of eireuhitiuj^ eiu^reiit obtained at that value of and Tth is the 
efi^eetive device temperature. Fittuig a nieaauretuejit of /\ Eia a fnnetiou of to E[[. ISIH 
combined with a Jmowledge of Ip, allows lis to extract Tn,. 

Meaaurejuent of Tn, was performed on two of the devices at each temperature setting. 
All average of at lea&t two independent lueiiHurements of the device temperature T^^ of each 
of two qubits is compared agauiat the mixing chamber thermometer temperatme reading 
(Tjv^Ar) hi Figure Uncertainty hi T,i^ was dominated by the uncertainty in the fit 

transition width for each measurement, which generally found to be larger than the 
stantlard (leviation of the separate uieasmementa. 
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time l^is) 

FIG. SS. The auuealiiig echedule is defined hy the applied fiux The qubits make a transition 

between heing uiouostable and bialable between toutrui points a = (10 jis, — 0.j92 and ^ = 
(158 ^9,-0.652 *a). 

Tho tojnpora,turn ostractHl Irotn MRT trannition rate widtlis {Tmur] is also plotted vh. 
Tfifxc fbr tomporatiurca below 40 mK, in Figure Ftom these plots it is dear that the two 
methodti generally a^;rc?e with eaeh other as wtiil aa with the :nisiug dianiber thermometer 
to within 3 inK over the temperatiire range used in the sxperiment. 

IV, QUANTUM SIMULATIONS 

To oljtJQii Letter quantitative underatandiiig of tlie beliavioiir of tlie system, a simulation 
waa eondutted to jnodel this experiment. The agreement between the numerical simulations 
can be aeen in panel (b) of Fig.|52] where both percentages [eKpcriment (theory)] are reported 
jiext to eat'h other for eaeh one of the low-energy eonformations. 

Our himnlation strategy is ay follows; We tirHt write a Hamiltonian for the fiupcrconduct- 
ing circuit based on standard circmt models for capaeitaiiees. uidnctanecSj and Joaephson 
junctiona. Tiiiti Hajniltonian is ejtpeeted to correctly dewcribe the behaviour of touplt^d 
rf-SQUIDs. We then numerically caJculato the evolution of the system based on tiiia Hjaiiil- 
tonian uning quantinn :nechaiueal etinations of motion which take into account couphng to 
an enviromueiit. Therefore, we predict the quantum evolutions for the same system Hjaiiil- 
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FIG. S^. Plots of (black circles) and T'.uiJT' [red triangles) vs. the temperature meB^ured with 
the Ruthemum Oxide thermometer mounted on the mixing diBjnber Tfuxc- 

touiaji, tlio saiiic! coupling Eo ciivirojiinciit, niid tht? saine type of noinc spi?ctral dt!Uhiti(?H. 
This provides a fair compEirisou to the cxperimcjital data. 



A. rf-SQUiD Haimltoniau 



Asiiuplificd \^rsioji of the rf-SQUID qubit used in our proct'SHor is illiJtitratod iu Fig. S5 
(A more complete desciiptioTi of the tictn-d qubita eaji be fouud iii Ref. It hsis two main, 
supcreondutting loops mid therefore two fltrx degrepK of frct>doin arid ^j. subject to 
external flux biases arid respectively. The Hamiltoiiian of such an rf-SQUlD is 
written ajs 

where Ci and Ca are parallel and series combinations of the junction capjicitaucea, gi and q-2 
are the sum and difference of the clmrgea stored in the two Josephaon Junctions respectively, 
and 

-2i:jCoa(7r*^/*u)cos{27f'&|/*,>), (S21) 
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iti A 2-dimt!riHioual pott?nljal with L, Lt^iii^ Uiu iiiductiuico^i pi the two loops aud $0— ''/^^t 
t}n: Hijx qiiiinliim. Wchavo assumed symmetric Joscphscia jmictions wfth Joscphsoo energies 
Ej= I^^{,/2iT, whurL' I^. lb tlio junctions' critical currc?nt. [A small asymmotry can be? tuii(?d 
away in nitii in the physical ijnplcnioiitatioiP.) 

At 0,)/2, the potential can bctonit? bistable ajid therefore form a two-dimensional 

double-well potential. If is sm-iJl cuough so that the deviation of '^^P Itoin be 
uegleotcd, then tbo two-dini(?Uhiouid chiasical poteutial U{^i,^^) can be approximated by 
a one-dimensio3ial double-well potential, -is shown in Fig. ^^ n. Howi.'ver, with our realistic 
qubit parameters, 0^ taunot be neglected and therefore Is accounted for in all om tiuinerical 
calculations. When = the two wells are symmetric with no energy bias between 

them. One can tilt the potential by ciianging Oi^, and eatabhbh an energy bifis, as depicted 
in Fi^. It is also possible to change the barrier height by changinj- 0^^. 

An array of such qubits can be modelled by auiuming contributionti of Eq. jjS20|| fro:n 
each device phis tenna that describe magnetic coupHng of the loops: 

Coupling between qubits i and j can be modelled ah a mutual inductance A/^j between loop 
1 of each pair of cuiipled qiibita: 



As discussed in Section 



lUC 



above, all paraiueteih, i.e., mdnctances La\ capacitances Cq\ 



and Josephson critical currents !c \ are measured independently for each qubit. 

To describe the system accurately we also need to introduce interaction with environment. 
Flux noise, which is the dominant noise in flux [[ubitn, couples to the ith qubit as fluctuations 
S^^l of the extennd flux a>L'i: 

Q=l 1 

The noise is much smaller foi the smaller loop than for the larger loop due to the 
loop siae. The flux noise S^aL is assumed to be uncorrelated between the qubits, which 
agrees with recent experimental observatioiP^. 
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1. Chip calibTatioji and device paTameter extraciioii 

Device pcir-ijiii^l.i.'fji were oxLr-icLod fui ihii rtiiinilalitmti ihruiigli i\ sicrit's of iadepcadent 
meaaurem^utH of qubit circiilEituig currt-ut, tuuiit?] Mplittiug A. find MKT ptaak spacing, A 
discussion of huw these uieasuremeuts are petfurriic?d is given in Rcf. |3i Parameter values 
used in siinulalions arc auiuniarised in T^bit^lTIIl 



B. Quantum Simuliition 

To siniulale lln.^ qufuituin LneciiiiiiitEii dynamics of the Hy&tem, we treat < ]S22[ |-( [524| as 
quantum meebanicai Hamiltoniaiis. Li that caaa, the charge is taken to be -m operator, 
whidi ia the momentum conjugate to the flux operator $1,'' with eommutatioii relatio:i: 
d\ ?q'] = iti. Unfortunately, It is impossible to calculate the dynamics of the system 
directly on the 2A-dimeuaion-d continuous poteuti^d quaiituui uieehajiieally. In&tead, we use 
energy diacretiKatioii as a meaiLs to simplify the caleulatio:i. The simplest way to accomplish 



this is to treat an rf-SQUID as a 2-atate syateni or qubit and replace ( S22 i by a coupled 
qubit Hamiltoniau. One may go further ajid keep more than two states per rf-SQlilD iji the 
calculation, as we shall discuss below. 



We first umuerically diagonalise the single rf-SQlJID Haiuiltoiiian [S20i to obtain the 
lowest Gigeiivaluea and eigenvectors. We treat the lowest few energy levels as the subapacc 
relevant for computation. We then write the Hamiltoniau in the basis of states that are 
localised within the wells. Such states are not true eigcjifunctiona of the Hamiltoniaji and, 
therefore, are metatitable towards tunnelling to the oppOHite well. HeneOj the remilting 
Harniltouian in such a basis will have olf-diagoiial tenns between states in the opposite 
wells but not between states withui eaeli well. The latter is because those states should be 
stationary^ within their own wells; any transition (relaxation) between them is only induced 
by the environment. 

Let \l) denote localised states withiTi the wells. We use even (odd) state numbers, i.e., 
I = 2ii (2jp-fl], with 11 = 0, 1; 2, .... to denote states that are locahsed in the left (right) 
well. For the lowest M energy levels (A/ is taken to be even], the cEcetive MxAl tmjuelhug 
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HEiEiiiltouiaji is written as 

Hs = J2^iM\^ A'a.>.^i(|2T0(2m+l| + |27n-fl){27i|) (S25) 

where; Et is the energy expectation value for state |^) and A'on.^nn-L in tlio tiinuelliug amplitude 
betwecTi states \2il) and |2m+l)| which exist in oppttaitc wells. Notice that there is no 
matrix element between states on the jiaine well; (2r;|//s|2m) = (2n+l|Hs|2m+l} = 0, 
which means that the states are metastable only towards tmmelling to the other aide, or 
the Htates are quasi-eigejistateH of the Hamiltojiian within their own sides. All parameters 
of the tumjelling Hainiltonian. i.e., Ef atid Kn- are extracted from the original rf-SQUlD 
Hamiltonian j ]S2[l[ | . For the 2-state qubit model we keep only the lowest two energy le^'els 



of [ S25 I. The eEeetive qubit Hamiltonian can be written as 



where 

We also go beyojid the 2-state model ajid keep 4 states per rf-SQlj'ID. Tliose 4 states can 
be represented by two coupled quhits, one of which represents the direction of persistent 
current or tlnx, and tiie other one generatuig intrawell ener^ levels, We represent the first 
(logical) qnbit b>^ Panli jnatricea Oa, and the extra (aneilla) qnbit by PauH matrices To- The 
effective HajuiltoJiian for those two coupled qubits can be written as 

It is CEisy to show that < ]S2S[ i ia equivalent^ up to a constant energy, to Hamiltoniaji j ]S25[ l . 
\\ith M = -L if 

€ = E^-Ei = E2- E-^. ojp = E.-Eu = E;h - E^, A = -2Koi, (329) 

= - A'ui ^ A'iHs ^■;J^ = 2A'fii = 2Ki2- (S30) 

As can be seen, the eoupli:ig bet\^^>n logical and aneilla qubits are of XX-hXZ type. Couiiling 
qubits to each other is accomplished using operators which represent the direction of the 
induced flux. The aneilla qubits remain uncoupled from each other and frojn other qubits. It 
slioiild be noted that the readout at the end of the evolution can only distinguiaii "left'" well 
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FIG. SiO. Parameters af the 4-level model for rf-SQlJID qubit ss a function of time during anneal- 
ing. 

from "right" wtiil hi tho doublu-vL^il potcutiEi] and cannot dLstinguish levels withui L'ath well. 
This is equivEilciit to reading out logic-d qubits and not ancilla qnbita, but as we mcntiont^d 
aboi^, only logical qublta c&iry information. 

To properly troat the ciiviromnoiit, wc oectl to write the intcraotiuu Haniiltoniiiji in 



the HubapacH o£ the lowest unergy levels iji terms of Pauli matrices. For quantum iiiiiulatioiiH 
we only consider noise coupling to the latgei loop in Fig. Let m consider a angle rf- 
SQUID and write the interaction Hamiltonian as 



(S31) 



Wt? define the? qiibit persistent current by 



I 



J, = -\{1\{^,-^^^)\1)\. 



(S32) 



Here, \vu take Ip to be independent of |/) for the low lying states cojisidcrcd, although in 
reality there could be a small dependence. The interaction H-uuiltonian can. then be written 



where 



l|Jf,nL|2) _ {l|fr,nt|3) 



{S33) 



{S34) 



The matrix elements (fl|i/ini|2) or (l|i/,ni,|3) are calculated directly via Eq. ^31 Uhiii^^ the 



dgenfunctioua of the rf-SQUlD thiniiltoniau. Eq. S20 ■. The values of Ip and X can therefore 
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be calculated Jiuiaericidly from the origiuEil rf-SQUlD Hajuilloiiiiui. Oiiiy Q [t'Liiajubi wliiti: 
shoiild bo characterised via its spectral doositvs which is the subject of Appcudis A. 

Quautum ewilution of the system was calculated iising a M-irkoviaii master equation for 
the dtiUHity matrix described iu Rcf. [Ti| Since the evolution i,s veiy slow {adiabaticj and 
temperature is low. only a aui-dl number of cjiergy levelti arc expected to be occupied diuiTig 
the ('\'olutiuri. We write the density juatrix in the instantaneous cnergj' eigenstate basis and 
trimcatc it tu the luwcst 24 euurgy levels, which waa foimtl to sufficiently dcacribe the type of 
evolution studied here. We utie both 2-Btate aiid 4-Btate models for rf-SQUlDs, aw described 
abova, in our simulations. The residt of the 4-level inoilel simulation is shown in Fig. [5^ . 



Appendix A: Noise spectral density 



The quantum noise operator Q = 2lpS^}r is related, to the flux noise as expected (for 
simplicity we only cojisider oue rf-SQlilD), ajid is clLaracterised by its correlation fmiction. 
Let us define the spectral density 

SH = r dt ,'"'{Q(t)im) = .i;;S*{^) (SI) 

J -X- 

where 

s*(u,) = r dt e'"{s^,,{t)s^,M) (S2) 

J—ca 

IB the spectral density of the flux noise. No direct measurement of 3n,{(J) at all frequencies 
IK available. We assume the spectral density is a sum of low frequency and high frequency 
comptiuenta: S-ti^) = SLF{i^) + ^hfCu^)- Ft>r the low frequeucy component we use 

{AyksT)ru^\..-\-^ 

with ci i^i 1, wliich at low uj behaves as 1/f noise: A^|[i;|~°. P-Lrameter A is measured from 
low frequeucy noise meiisuremenrf^ and is found to he n'^^,. 

The liigh frequency parts of the spectral density is assuuied to be ohmic, 

where is the upper cutoff frequency, ij is the dimensiojiless coupling coefficient, and /po 
is the value of persistent cmrent at which tj is measured. The couphng coefficient and the 
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persistent mrrt^ut iirc luiiiid, uaitig Matrostopic roaouuiit tiiiiucllitig oxpt^riuioiit (MKT), to 
be ij ^ 0.1 ajid /p,, ^ 1 fiA. The dctailH of extraction of via MRT arc preacntcd elaewhcrd^. 
This IciivoH iiu free pariuiieterH for the (juaiituiLi siuiulatiuiiH. 
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